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X 



Context. Blazars are radio-loud active galactic nuclei (AGN) with a jet pointing at small angles towards the observer The overall 
emitted spectrum is typically non-thermal, and in some cases the emission and/or absorption lines are so faint as to prevent the 
determination of the redshift based on optical spectroscopy methods. PKS 0447-439 is a bright blazar recently detected at very high 
energy. The redshift of the source is still disputed: a recent spectral analysis reports only a lower limit of z > 1.246, which contradicts 
the previous measure of z = 0.205 reported in the literature. 

Aims. We aim to give a redshift estimate of the blazar PKS 0447-439 based on combined GeV (Fermi/LAT) and TeV (H.E.S.S.) 
observations. 

Methods. Taking into account the absorption of TeV photons by the interaction with the extragalactic background light (EBL), we 
applied the method we developed in a previous work to derive the redshift of PKS 0447-439. Moreover, we compiled the overall 
spectral energy distribution (SED) using optical-UV, soft X-ray, and y-ray data, nearly simultaneous to the H.E.S.S. observations 
at TeV energies. Finally we modelled the spectral energy distribution (SED) within the framework of a homogeneous, leptonic 
synchrotron self-Compton (SSC) model. 

Results. Using the recent TeV spectrum measured by H.E.S.S. we obtain for PKS 0447-439 a redshift of Zrec = 0.20 ± 0.05, which is 
our estimate on the source distance. This value agrees very well with the value reported in the literature and confirms that our method 
can be successfully used to constrain blazars distances. Assuming this distance, the SED can be well fitted with the above mentioned 
model. The physical parameters that we find suggest a strongly matter-dominated jet. 

Conclusions. Our analysis confirms that the redshift of PKS 0447-439 is likely 0.2, and supports the result present in the literature. 
Key words. BL Lacertae objects: individual: PKS 0447-439 - Galaxies: distances and redshifts - Gamma rays: galaxies 



1. Introduction 

The number of known AGN whose emission extends up to the 
very high energy (VHE; E > 100 GeV) y-ray band has more 
than doubled in the last couple of years. Nowadays, the TeV 
sources catalogue counts 44 located both in the north- 

ern and in the southern hemisphere. This achievement has been 
made possible thanks to the good performances of the last gen- 
eration of imaging atmospheric Cherenkov telescopes (lACT), 
namely MAGIC, H.E.S.S. and VERITAS. A key ingredient for 
the detection of new sources has been the cooperation of these 
telescopes with satellite experiments, especially those operating 
at optical. X-ray and soft y-ray frequencies. 

The great majority of AGN detected at VHE belongs to the 
class of blazars. In these sources, the emission is produced in 
relativistic jets that point towards the observer. 

PKS 0447-439 is a blazar located at RA(J2000) = 
04^9^24^88, Dec.(J2000) = -43°50'09'.'7. It was discov- 
ered in 1981 at rad io wavelenghts by the Molonglo Telescope 
(iLarge et all Il98lh and de tected by the PMN radio Survey in 
1993 (iGreeorv et aIi[T99l . 

The source was sub sequently detected in the U V with EUVE 
(iLampton et all Il997h and X-rays with ROSAT dWhite et all 



See TeVCat (http://tevcat.uchicago.edu/ 1 and 
(http://www.mpp.mpg.de/~rwagner/sources/) pages. 
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[1994 1). Some confusion arose in the determination of the opti- 
cal counterpart because at first a near-UV bright spectrum with 
prominent emission lines in the optical was related to PKS0447- 
439 (Craig et al., 1997 ). The source was therefore c l assifie d as 
a Type I Seyfert at z = 0.107. Instead iPerlman et all (Il998h re- 
ported a featureless optical continuum typical of a BL Lac ob- 
ject. 

We investigated this classification discrepancy, also inspect- 
ing some archival Swift/WOT exposures of the field. We un- 
derstood that the association of the Seyfert spectrum with this 
source is probably due to some databasing mistake, because the 
PKS 0447-439 coordinates reported in Craig et al.1 (Il997h are 
incorrect, and point to a region that is source-free both in their 
optical finding chart and in our UVOT images. 

The reds h ift of PKS 0447-439 was measured by 
iPerlman et al.l (Il998h . who reported a value of z ~ 0.205 
based on few weak absorption features in the optical spec- 
trum that were interpreted as th e Cali,K doublet. Another study 
performed by iLandt & Bignalll (12008 ) provided only a lower 
limit of 0.176, based on the non-detection of the host galaxy, 
whi ch is assumed to be a g iant elliptical of inferred luminos- 
ity (iPiranomo nte et all. '2007). Recently, a different lower limit 
z > 1.246 was obtained (Landt, 2012) based on the analysis of 
absorption lines in the optical spect rum. This new value is well 
above the measurement reported bv lPerlman et al.l (Il998l) . 
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The Large Area Telescope (LAT) instrument on board Fermi 
identified the source as one of the brightest of the southern 
hemis phere in the 100MeV-300GeV energy range (Abdo et al., 
I2009h . Dedicated studies on the weekly light curve above 
300 MeV revealed th at the source is variable in this energy range 
dAbdo et al.L 1201 Oal) . Inspired by the FermifLA^ observations, 
the H.E.S.S. lACT observatory performed a dense observation 
campaign between November 2009 and January 2010. This led 
to the detection of a VHE signal, statistically significant above 
the 13 cr level and position ally consistent wit h PKS 0447-439, 
as preliminarily reported in lZech et al.l (1201 ll) . 

In this Paper, we derive an independent estimate of the 
distance to PKS 0447-439 using combined GeV and TeV 
spectral information , following the method we proposed in 
iPrandini et al.l ( l2010l) . Our goal is to provide a new measure- 
ment of the source's redshift, which is still uncertain, as dis- 
cussed above. In the next section we briefly outline the method 
and present the r esult, which basically confirms the estimate of 
z ~ 0.2 given bv iPerlnian et al.l (119981) . Then, we build a quasi- 
simultaneous SED, which is fitted with a homogeneous, leptonic 
SSC model. This allows the determination of the main physical 
parameters governing the jet physics. 

2. Inferring the distance of TeV emitting blazars 

The VHE spectrum of blazars suffers from the absorption aris- 
ing fro m the interaction with the extragalactic background light 
(EBL) (Hau ser & Dweld 1200 Ih . This absorption, caused by 
electron-positron pair production (jNikishov, 1962), induces a 
partial/total deformation of the VHE part of the spectrum that 
is strongly redshift-dependent. In general, for nearby sources 
that are located at redshift below 0.1, it affects the spectral 
points above some TeV. At higher redshifts, between 0.1 to 
0.5, it affects the spectrum already at some hundred of GeV, 
while above 0.5 it becomes effective already above some GeV. 
Therefore, given an intrinsic VHE spectrum, the observed one 
depends on the distance of the emitter In other words, the 
optical depth associated to the absorption process depends on 
the distance covered by the energetic photon and on its en- 
ergy. Unfortunately, owing to strong foreground emissions that 
are difficult to suppress, there are no solid measurements of 
EBL. Upper and lower limits are provided by indirect tech- 
niques such as galaxy counts, which provide solid lower limits. 
In addition to these measurements, many models of the EBL en- 
ergy density and evolution have been proposed in the last years 
(IStecker et al.'.'2006^.'Franc eschini et all 120081: iKneiske & Dolel 
[2OIO; Dommauez et al., 20 lH). In this work we adopt the model 
presented in Francesc hini et al.l (12008'). 

To estimate the distance to PKS 0447-439, w e applied the 
m ethod that we propos ed in iPrandini et al.l (l2010h and updated 
in iPrandini et al.l (1201 ll) . This method is based on the compari- 
son between the high energy (HE; 0.1 < E < 100 GeV) y-ray 
spectrum measured by Fermij'LP^I and that measured at VHE 
by the last generation of imaging atmospheric Cherenkov tele- 
scopes (MAGIC, H.E.S.S. and VERITAS). 

W e refer the reader to the original paper (IPrandini et al.L 
I2OIOI) for the full description of the method. Briefly, to infer 
the distance of an HE and VHE gamma-ray emitting blazar, we 
calculate the redshift, z*, at which the power law slope fitting 
the VHE spectrum corrected for the EBL absorption equals the 
slope measured by Fenni/LAJL at lower energies. As empirically 
demonstrated in the paper using a set of known redshift sources, 
this value is related to the true redshift of the source by a simple 
linear relation. Therefore, once one has obtained z*, it is possi- 
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Fig. 1. Observed (light grey, open circles) and deabsorbed VHE 
spectra of PKS 0447-439, as report ed in Zech et al. ( 201i]) 
Deabsorption is applied assuming the Franceschini et al.l (l2008l) 
EBL model, and a redshift of the source z* - 0.36 (black, filled 
circles) and Zrec = 0.2 (dark grey, filled squares). 



ble to give an estimate on the source distance, Zmc, by using the 
inverse formula: 



with A and B given in IPrandini et al.l (1201 ll) . The law was 
applied to infer the distance of the unknown redshift blazar 
PKS 1424-1-240. The value obtained has been confirmed in a 
study on the host galaxy performed bv Meisner (.2010.) . 

2. 1. The reconstructed z of PKS 0447-439 

The blazai- PKS 0447-439 was observed by the H.E.S.S. lACT 
observatory between November 2009 to January 2010, for a total 
of 13.5 hours of data. A clear signal of 13.8 cr of signifi cance 
was detected at energies above 250 GeV dZech et 3111201 ll ). The 
preliminary spectrum extends from 300 GeV up to more than 
1 TeV, and is compatible with a steep simple power law of index 
F = 4.36 + 0.49, where the error reported is statistical only. The 
significant spectral points are drawn in Fig.[r| open markers. 

The Fermi/LAT spectral slope in the energy range 0.1- 
lOOG eV is Flat ^1.95 ± 0.03, as reported in the IFGL cata- 
logue dAbdo et al.L 1201 Obi) . We adopted this slope, and not the 
one extracted from only simultaneous data, to be consistent with 
the analysis performed for the determination of the z* - Zme law. 

The redshift that we obtain when requiring that the deab- 
sorbed PKS 0447-439 spectrum has a slope equal to Flat is 
z* =0.357 + 0.065, where the error takes into account both the 
errors on the TeV slope and on the FermifLP^ slope. The cor- 
responding spectrum is drawn in Figure [Tj filled black circles. 
A simple power law fits the data very well, with a probability 
higher than 99%, and ^^/ndf = 0. 1/3. 

F igure [2] shows t he z* - Zmie plot in linear scale, adapted 
from IPrandini et aP (l201ll) . The grey markers represent the z* 
values of the sixteen sources with known redshifts used to infer 
the linear relation (dashed-dotted line). The continuous line is 
the z* = Zinie line. All points lie above or on the line, meaning 
that the values z* are higher than the real redshifts. This indi- 
cates z* as a good upper limit on the source distance. In our case, 
we have z* of 0.357 + 0.065 and the corresponding upper limit 
is 0.49 at 2 sigma level. Applying Equation (Tj to our data, we 
obtain for PKS 0447-439 the distance z„, = 0.20 + 0.05, where 
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Fig. 2. Diagram adapted from iPrandini et al.l (1201 ih represent- 
ing the linear relation between the redshifts of known dis- 
tance sources and their z* values. The value of z* obtained for 
PKS 0447-439 in this work is superimposed in black. 

the error is statistical only. The differential energy spectrum ob- 
tained assuming a distance z of 0.20 is drawn in Figure [T] filled 
squares. It lies between the observed spectrum (not corrected for 
EBL absorption) and the spectrum obtained assuming a distance 
z*. A simple power law of index F = 3.2 + 0.5 fits the data very 
well (probability > 99%). 

The estimated distance obtained perfectly a grees with the 
value z = 0.205 reported bv lPerlman et alJ (Il998h . and also con- 
firms the lower limit of z > 0.176 based on photometric es- 
timates. Therefore, our study strongly supports the value of 
z = 0.205 for PKS 0447-439. 

3. Spectral energy distribution 

The single-epoch multiwavelength SED of PKS 0447-439 is 
plotted in Fig. |3] and includes, in addition to the H.E.S.S. 
data presented in previous section, Swift/UVOT, Swift/XRT, and 
FennifLALT data acquired during the epoch of the H.E.S.S. de- 
tection, and other historical data. 

Swift observed PKS 0447 -439 for about o ne week at the 
end of t he H.E.S.S. campaign (IZech et al.L 1201 Ih . Data from the 
UVOT (Roming et al., 2005) observation taken on December 
25 (obs. ID 00038100009) were analysed by means of the 
uvotimsum and uvotsource tasks with a source region of 5", 
while the background was extracted from a concentric source- 
free annular region with inner and outer radii of 10" and 15", re- 
spectively. The extracted magnitudes were corrected for Galactic 
extinction using th e values of Schlegel et al. (.1998.) and apply- 
ing the formulae bv lPeil (Il992h for th e UV filters, and ev entually 
were converted into fluxes followingl Poole et"al] (|2008|) . 

XRT (Burrows et al., 2005) vFy points corresponding to the 
observation of Dec 25 (MJD 55190, obs. ID 00038100010) 
were obtained thro ugh the ASl/ASDC on-line analysis tooQ 
(IStratta et al.L l2010l) . The X-ray fluxes were obtained assuming 
a power law spectrum (best-fit photon index F - -2.46) and 
fixing the absorption to the Galactic value A^^"' - 1.24 x 10^" 
cm-2 (after lKalberla et al.Ll2005b . 

We retrieved the publicly available data tak en by the LAT y- 
ray telescope on board the Fermi satellite ( At wood et all l2009t) 
in scanning mode from the NASA databas^. We selected the 

^ http://tools.asdc.asi.it/ 
^ http://fermi.gsfc.nasa.gov/ 



good quality ("DIFFUSE" class) events observed within 10° 
from the source position, taken between MJD 55170 and MJD 
55190, and with measured energy in the 0.2-lOOGeV interval. 
We excluded events observed at zenith distances greater than 
105° to avoid contamination from the Earth albedo. We per- 
formed the analysis by means of the standard science tools, 
V. 9.18.6, including Galactic and isotropic extragalactic back- 
grounds and the P6 V3 DIFFUSE instrumental response func- 
tion. We applied an unbinned likelihood algorithm (gtlike) to 
the data, modelling the source spectrum with a power law model, 
with the integral flux in the 0.2-lOOGeV energy band and pho- 
ton index left as free parame ters. W e clearly detected the source 
with a test statistics (Matt ox et al.L 11996) TS = 154, and de- 
termined Fo.2-iooGeV = 6.93 + 1.26 X 10"*^ phcm^^s"' and 
F = 1.93 + 0.14. 

We also performed the same analysis on the dataset split into 
logarithmically equal bins of energy, accepting spectral points 
with TS > 25 and more than three events attributed to the source 
in the model. The results are plotted in Figure|3] 

PKS 0447-439 belongs to the group of highly peaked BL 
Lac objects. Within the framework of leptonic models, the emis- 
sion of this type of sou rces is generally modelled with the one- 
zone SSC model (e.g. Tavecchio et al 1 119981 l2010HAbdo et al.L 
2011). The idea that the emission mainly originates in a sin- 
gle, uniform, region is supported by the correlated variability 
observed at different frequenc ies, especially X-rays and gamma- 
rays (e.g. lFossatietal.Lll998h . 

The good spectral coverage allows one to firmly constrain 
the low-energy bump, which is related to synchrotron radia- 
tion produced by accelerated electrons in the jet, whose peak 
lies around 10"' Hz (see Fig. 3). LAT and H.E.S.S. data track 
the so-called SSC component quite we ll, which shows a broa d 
peak in the GeV region. As detailed in lTavecchio et al.l ( Il998l) . 
in the framework of the one-zone SSC model knowing of the 
peak frequencies and luminosities, together with an estimate of 
the variability timescale, allows one to determine the physical 
parameters of the emitting region. An indication of the variabil- 
ity timescale can be derived from the multi-frequency observa- 
tions, especially those at X-ray energies, which reveal variability 
on timescale of days. We therefore assume an upper limit of 
fvai- < 1 day for the variability timescale. 

We modelled t he SED w ith the one-zone leptonic model de- 
scribed in Maraschi & Tave cchio (2003). The emission region is 
spherical with radius R, in motion with bulk Lorentz factor F at 
an angle 6 with respect to the line of sight. Special relativistic 
effects are fully accounted for by the relativistic Doppler fac- 
tor, 6 - [F(l -jScos0)]"'. The energy distribution of the emit- 
ting electrons is assumed to be well described by a smoothly 
connected broken power law function, with minimum, maxi- 
mum and break Lorentz factor jmin, 7max and yb, respectively. 
The SSC emis sion is calcula ted assuming the full Klein-Nishina 
cross section (I.TonesL Il968l) . We recall that within the frame- 
work of the one-zone SSC model the synchrotron self-absorption 
causes the source to be opaque below frequencies of about 
10" - 10'^ Hz. The emission below these frequencies is there- 
fore produced by more distant, transparent regions of the jet. 

Two possible SEDs are reported in Fig. 3, while the corre- 
sponding input parameters are listed in Table [T] In the table we 
also report the derived powers carried by the different compo- 
nents, relativistic electrons, magnetic field and protons (assum- 
ing a composition of one cold proton per relativistic electron) 
and the total radiative power of the jet. Pi ^ Lobs/^^ (in which 
we assumed 6 ~ F). 
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Log u [Hz] 

Fig. 3. SED of PKS 0447-439 during the epoch of H.E.S.S. observation (November 2009 - January 2010). Red points represent 
Swift/UWOT, Swift/XRT, FermijhP^ and de-absorbed (assuming z - 0.2) H.E.S.S. data. Cyan data are the observed H.E.S.S. data. 
Green symbols report historical data obtained through the ASl/ASDC tools. Green LAT points are from the ILAC catalogue. The 
black lines are the result of the one-zone emission leptonic model discussed in the text for two different values of the minimum 
Lorentz factor of the emitting electrons, y^in = 1 (solid) and ymin = 3 x lO^* (dashed). See text for details. 
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Table 1. Input model parameters for the models reported in Fig. 3. We report the minimum, break and maximum Lorentz factor and 
the low- and high-energy slope of the electron energy distribution, the magnetic field intensity, the electron density, the radius of 
the emitting region and its Doppler factor. We also report the derived power carried by electrons, magnetic field, protons (assuming 
one cold proton per emitting relativistic electron) and the total radiative luminosity. 



The jet power is strongly dependent on the total density of 
particles in the jet that, in turn, is dominated by the number of 
relativistic electrons at jmrn- While this parameter can be rela- 
tively well determined for FSRQs, because the low-energy end 
of the electron energy distribution is directly accessible through 
the inverse Compton emission that falls in the X-ray band (e.g. 
iGhiselUni et aliuOlOh . for BL Lac objects emitting through SSC 
the constraints are quite loose. As an example of the impact 
of different values of y^in on the derived SED, we report in 
Fig. [3] two curves, corresponding to y^in = 1 (solid line) and 



Train = 3x10^ (dashed). The two curves clearly show that the 
value of yniin mainly affects the low-energy branch of the syn- 
chrotron and SSC bumps. Clearly, any value between 1 and 10^ 
is allowed by the available data, which determines a large uncer- 
tainty on the inferred powers (see Table 1). 

The magnetic field energy density appears to be several or- 
ders of magnitude below the equipartition value with the electron 
energy density and, consequently, the corresponding Poynting 
flux is negligible compared to the power carried by protons 
and electrons. This is a feature that PKS 0447-439 shares 
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with several other Te V BL Lacs (e.g. iTavecchio et al.L l2010t 
iGhisellini etani2010h . 

4. Conclusions 

We have presented the results of our study on the distance 
to PKS 0447-439. This work was triggered by the prelimi- 
nary measurement of the differential energy spectrum emitted 
by the source at VHE gamma-ra ys performed by the H.E.S.S. 
Collaboration dZech et al.L 1201 lb . The redshift we found with 
our method, based on both HE and VHE s pectra, perfectl y agree s 
with the measurement performed by iPerlman et al.l (Il998h . 
which reports a redshift of 0.205. But it contradicts with the re- 
cent assertion of Landt (2012), who claimed a redshift higher 
than 1.246. 

Assuming this redshift, we have built a broadband SED in- 
cluding quasi-simultaneous ultraviolet. X-ray and gamma-ray 
data acquired during the epoch of the H.E.S.S. detection. The 
SED can be well fitted within the framework of a homogeneous, 
leptonic SSC model. The available data leave a large uncer- 
tainty on one of the parameters of the model, jmin, the minimum 
Lorentz factor of the emitting electrons. This uncertainty propa- 
gates to the inferred powers. In any case, the physical parameters 
that we found suggest a strongly matter-dominated jet, in agree- 
ment with the other sources observed at these extreme energies. 

The method used in this paper has great potential. A large 
part of GeV emitting blazars, indeed, about 60% of the BL Lac 
objects prese nt in the second catalogue of AGN detected by 
Fermi/LACr dAckermann et al.L [201 Z) . has unknown redshift. If 
detected in the VHE band, the redshift of these powerful objects 
could be investigated through this technique. Therefore, we en- 
courage the observation of promising TeV emitters among the 
unknown/uncertain redshift Fermi blazars. 
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